Multi-Orbital Quantum Phase Diffusion 
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The collapses and revivals of a coherent matter wave field of interacting particles can serve as 
a sensitive interferometric probe of the interactions and the number statistics of the underlying 
quantum field. Here we show how the ability to observe long time traces of collapse and revival 
dynamics of a Bose-Einstein condensate loaded into a three-dimensional (3D) optical lattice allowed 
us to directly reveal the atom number statistics and the presence of effective coherent multi-particle 
interactions in a lattice. The multi-particle interactions are generated via virtual transitions to 
higher lattice orbitals and can find use for simulations of effective field theories with ultracold atoms 
in optical lattices [Ij . We measured their absolute strengths up to the case of six-particle interactions 
and compare our findings with theory. 



Coherent quantum states represent the most robust 
and stable field solutions in physics [2]. Character- 
ized by a single amplitude and phase, they have found 
widespread use in the quantum description of classical co- 
herent fields, ranging from laser light to coherent matter 
waves in superconductors, superfluids or atomic Bose- 
Einstein condensates. Whenever interactions between 
the underlying particles are present, or - more generally - 
whenever the phase evolution of the number states (Fock 
states) that form the coherent state is nonlinear in par- 
ticle number, coherent states can undergo an intriguing 
sequence of collapses and revivals. In such a sequence, 
the quantum state first evolves into a highly correlated 
and entangled state where at the time of the collapse the 
classical field vanishes, however, at a later time the en- 
tanglement is unraveled again and the original classical 
field is ideally recreated. 

Remarkable examples of such collapses and revivals 
have been observed for a coherent light field interacting 
with a single atom in Cavity Quantum Electrodynamics 
H 0], for a classical oscillation of a single ion held in a 
trap [5J [B] or for a matter wave field of a Bose-Einstein 
condensate [7HHH]. In the latter case, where the dynam- 
ics are induced by non-linear two-particle interactions 
between the atoms, the underlying process of quantum 
phase diffusion [H El El El El ES] has been proposed 
to be useful for the creation of highly correlated and en- 
tangled quantum states, such as large Schrodinger cat or 
spin-squeezed states [HEZlIllli the latter of which have 
been demonstrated recently [19]. In all these cases, it is 
typically assumed that the atoms occupy a single spa- 
tial orbital of the system. Atom-atom interactions can 
however promote particles to higher lying orbitals. The 
admixture of these excited orbitals modifies the shape of 
the spatial wavefunction and gives rise to renormalized 
interaction energies depending on the atom number, in 
close analogy to the so-called configuration interaction in 



quantum chemistry [20 . An instructive way to envisage 
the physical processes is provided by effective field the- 
ory: Virtual transitions of atoms to excited orbitals gen- 
erate effective multi-particle interactions as higher order 
corrections to the single-orbital two-particle interaction 
[I] . Previous experiments have only allowed the observa- 
tion of few cycles of quantum phase diffusion dynamics 
d H [H| an d therefore could not reveal the striking signa- 
tures and consequences of multi-orbital effects, which are 
expected to lead to a multitude of novel many-body quan- 
tum phases [2H 122 • Here, we have been able to observe 
more than 40 collapses and revivals in long time traces 
of the quantum evolution of coherent matter wave fields 
trapped in a 3D optical lattice potential. Multi-orbital 
effects give rise to a discrete set of distinct two-particle 
interaction energies for different Fock states constitut- 
ing a coherent field. These in turn result in character- 
istic beat signals in the collapse and revival time traces, 
from which we infer the occupation of Fock states, mea- 
sure their energies with high precision and determine the 
strength of the effective multi-particle interactions. Re- 
cently, multi-particle interactions have been identified via 
loss resonances for three- and four-body recombination 
|23ll24ll2"5] . however coherent multi-particle interactions, 
which we demonstrate here in the form of effective inter- 
actions, have so far defied observation. 

Let us consider a single site of a deep optical lattice 
filled with n atoms occupying the lowest single-particle 
state of the system V'o( r )- Assuming weak interactions 
and excluding multi-orbital effects, collisions between the 
atoms lead to an energy shift of the single-orbital Fock 
state energy given by = Un{n — l)/2 represent- 

ing the onsite interaction energy of the Bose-Hubbard 
model. Here U — 4Trh 2 a s /m J |^o(r)| 4 d 3 r denotes the 
two-particle interaction energy, determined by the onsite 
wavefunction ipa(r), the s-wave scattering length a s and 
the mass of an atom m. Within the restriction to the 
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FIG. 1: (A) A Bose-Einstein condensate loaded into a weak 
optical lattice forms a superfluid with each atom being de- 
localized over several lattice sites. The quantum state on 
each site can be expressed as a superposition of Fock states 
\n) with amplitudes c n . (B) For repulsive onsite interac- 
tions virtual transitions to higher lattice orbitals broaden the 
ground state wavefunction on a lattice site compared to the 
non-interacting system (gray dashed line). This gives rise 
to characteristic Fock state energies, which can be described 
by effective multi-particle interactions. (C) A coherent state 
of interacting atoms confined to a deep lattice well under- 
goes multi-orbital quantum phase diffusion (blue solid line). 
The beat signal indicates coherent multi-particle interactions. 
The dynamics are markedly different from the monochromatic 
evolution expected in a single-orbital model with a single two- 
particle interaction energy U (gray solid line). 



lowest vibrational state of the system, U is independent 
of the filling n at the lattice site. 

Inducing virtual transitions to higher vibrational 
states, interactions modify the shape of the ground state 
wavefunction (Fig. |TJ3) and the two-particle interac- 
tion energy U itself becomes atom number dependent 

[TlllDllBHTlliEllSH]- In this case ' the multi -° rbital Fock 

state energies in a single lattice well can be approximated 
by the expansion 



E%° = - U 2 n(n - 1) + - U 3 n(n - l)(n - 2) + 
t^4n(n-l)(n-2)(n-3) + ..., 



(1) 



representing the eigenenergies of the corresponding effec- 
tive Hamiltonian [1]. Here, coherent multi-particle inter- 
actions become explicitly visible with the characteristic 
strength of the m-particle interactions being denoted by 

Urn • 

An efficient way to experimentally probe the eigenener- 
gies of a Hamiltonian, is to monitor the non-equilibrium 
dynamics of a quantum state prepared in a superposition 
of different eigenstates. In our case, the state is a super- 
position of different atomic Fock states \n) and can be ex- 
pressed in the general form \i/j(t)) = X^^Lo c n e~ lEnt ^ h \n) . 
Experimentally, a 3D array of such states is created 
by loading a BEC into a shallow 3D lattice potential 
(Fig. [l]^-). The time evolution of the ensemble can be 
probed by analyzing the visibility of the atomic inter- 
ference pattern as observed after rapid switch-off of the 
lattice potential and subsequent time-of-ffight expansion 
[30] . For an array of identical states \tp(t)), the visibil- 
ity of the interference pattern relates monotonically to 
|(f/'(t)|a|V'(i))| 2 , where a denotes the annihilation oper- 
ator at a lattice site. Hence, the dynamical evolution 
of the matter wave field, the quantum phase diffusion of 
\ip(t)), is given by: 



\w)\*w))\ 2 = m\ 2 = 



oo 

E 
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To illustrate the essential dynamics of Eq. ([2]), we as- 
sume the initial state \tf>) to be a coherent state with 
c n = e~l Q l I 2 where a — Jn ■ e 1 ^ denotes the com- 
plex field amplitude with the mean atom number n and 
initial phase <fr. In a single-orbital model with eigenener- 
gies E%° = Un(n — l)/2, periodic collapses and revivals 
of the visibility are expected, which can analytically be 
expressed by |(a)| 2 /n = e 2n ( cos {Ut/K)-i) featuring a single 
frequency U/h and its higher harmonics (Fig. [TfU (gray 
solid line)) [H[T5]. 

Multi-orbital effects, however, reach beyond the pic- 
ture of monochromatic matter wave dynamics: the time- 
evolution of | (a) | 2 contains multiple frequency compo- 
nents since Fock state energies are no longer integer mul- 
tiples of U, as shown in Fig. [Tp (blue solid line). As- 
suming repulsive interactions, we expect frequencies of 
order 0(U/h) to be given by E 2 /h > (E 3 - 2E 2 )/h > 
(Ei - 2E 3 + E 2 )/h > (E 5 - 2£ 4 + E 3 )/h > ... as well 
as higher order frequency components 0(t ■ U/h), with I 
being a positive integer. 

The detection of multi-orbital quantum phase diffusion 
over long times results in a high spectral resolution, al- 
lowing for a precise measurement of individual Fock state 
energies E^f° and multi-particle interaction strengths. 
In addition, Eq. ^ implies that the spectral weight of 
each contributing frequency is directly connected to the 
coefficients c„ of the initial state Hence, the sta- 

tistical contribution of single Fock states to the global 
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FIG. 2: Multi-orbital quantum phase diffusion of an atom number superposition state. (A) The collapse and revival 
dynamics of a number-squeezed superposition state in a deep optical lattice is shown. A BEC of about 1.9 ±0.3 x 10 Rb atoms 
was adiabatically loaded to a Vl — 8 E rec lattice within 100 ms. Phase evolution was induced by a non-adiabatic jump into a 
Vh = 41 E rec deep lattice, preserving superposition states with finite number fluctuations and an ensemble averaged (central) 
mean atom number of about (n) = 1.0 (n c = 2.5). Simultaneously with the lattice jump, the underlying harmonic confinement 
was instantaneously minimized (C) to optimize the coherence time. The phase evolution shows a beat note signature resulting 
from coherent multi-particle interactions with different interaction strengths. Each data point corresponds to a single run of 
the experiment. The solid line interpolates the data and serves as a guide to the eye. (B) Spectral analysis of the time trace 
(A) using a numerical Fourier transform reveals the contributing frequencies. The solid line shows Gaussian fits to the peaks. 



superfluid quantum many-body state can be obtained 
complementary to previous work [3"T1 I3"2l |3"3"] . 

Our experiments began with an atomic Bose-Einstein 
condensate of repulsively interacting 87 Rb atoms in the 
\F = l,m,F = +1) hyperfine state, with variable atom 
numbers between 1.2 x 10 5 and 3.5 x 10 5 . The atoms 
were initially held in a pancake-shaped crossed optical 
dipole trap at a wavelength of Xdi P — 1030 nm with trap 
frequencies tu z = 2tt x 130 Hz in the direction of gravity 
and uj± — 2tt x 32 Hz in the orthogonal plane. Subse- 
quently a 3D optical lattice (A = 738 nm) of simple cubic 
type was adiabatically ramped up to lattice depths Vl 
between 3 — 13 E rec , where E rec = h 2 /(2m\ 2 ) denotes 
the recoil energy. For these lattice depths the many- 
body ground state of the system is a superfluid and we 
expect the states on a lattice site to range from coher- 
ent states for shallow lattice depths to highly number 
squeezed states for deeper lattices [7j [1J3 [3TJ [33] . 

A sudden increase of the lattice depth from Vl to Vjj 
ranging between 25—41 E rec strongly suppressed the tun- 
nel coupling and essentially froze out the atom number 



distribution on each site. In this regime, the time evolu- 
tion of each site is governed by the Fock state energies of 
Eq. (JTJ and the quantum phase diffusion process started. 
The jump to Vh was carried out within 50 /is which 
is slow enough to avoid populating higher vibrational 
states, but fast compared to tunneling within the first 
band. Previous experiments were limited to follow quan- 
tum phase diffusion only for short traces with a maximum 
of 4-5 collapse and revival cycles [7J [51 H] since a global 
harmonic confinement led to rapid relative dephasing of 
lattice sites. Here, the different detunings of the optical 
lattice (blue-detuned) and the dipole trap (red-detuned) 
with respect to the atomic resonance (A res = 780 nm) al- 
lowed us to change the underlying harmonic confinement 
during the experimental sequence (Fig. [2p). Along with 
the jump to Vh we instantaneously reduced the dipole 
trap in order to cancel harmonic confinement in the hor- 
izontal plane, boosting the coherence time of quantum 
phase diffusion. Under such conditions, the atomic sam- 
ple stayed trapped due to the presence of the deep opti- 
cal lattice Vh- After letting the system evolve for a hold 
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FIG. 3: Multi-orbital energies and multi-particle interactions. (A) Long collapse and revival traces were recorded 
under identical loading conditions (Vl = 8 E rec ), but variable lattice depths Vh during phase evolution. Numerical Fourier 
transform of the time traces reveals the contributing frequencies of orders 0(U/h) (red circles) and 0(2U/h) (blue circles) 
with a typical uncertainty of ±50 Hz, the shading of the data points reflects the relative spectral weight. The solid lines with 
gray shading denote the theoretically expected frequencies at an s-wave scattering length a s — 102 ± 2 ao as derived for a 
basis set with 4 3 orbitals. Calculation on a smaller basis set with 3 3 orbitals yields slightly higher energies (dashed lines). 
The black dotted lines show the single-orbital interaction energy of the Bose-Hubbard model (BHM). At low lattice depths, 
only the strongest peaks can be resolved mainly due to smaller peak spacings. (B) Effective two-, three- and four-particle 
interaction strengths as derived from experiment and theory (a s = 102 ±2 ao, 4 3 orbitals). The error bars correspond to 
standard propagation of uncertainty. 



time t, all trapping potentials were switched-off simulta- 
neously and an absorption image of the matter wave in- 
terference pattern was recorded after 10 ms time-of- flight. 
The visibility of the interference pattern was evaluated, 
being a measure for the ensemble average of |(a)| 2 [3"U] . 

A typical time trace (Vl = 8 E rec ) of quantum phase 
diffusion is shown in Fig.[2jY displaying about 40 revivals. 
On top of a fast series of collapses and revivals, we ob- 
serve a slower modulation of the envelope indicating a 
beat between at least two different energies in the sys- 
tem. The spectral content of the time trace has been 
quantitatively analyzed using a numerical Fourier trans- 
form of the data (Fig. |2j3). We clearly identify three dis- 
tinct frequency components of order 0(U/h) present in 
the time trace, the smallest one originating from sites 
occupied by up to four atoms. It is remarkable that 
our technique can reveal even very small Fock state am- 
plitudes c n due to a heterodyne effect with other Fock 
states |n — 1) and \n — 2) of typically larger amplitudes 
c„_i and c„_2 (Eq.[2]). The residual damping on the time 
trace most likely results from residual harmonic confine- 
ment (along the z-direction) , residual tunneling [34] via 
higher bands and the finite extension of the atomic en- 
semble, sampling a slightly inhomogeneous distribution 
of lattice depths due to the Gaussian shape of the lattice 
laser beams. 

In order to map out the dependence of the multi-orbital 
Fock state energies E^ 4 ° on lattice depth, we recorded 
several long collapse and revival traces with Vh between 



25 — 41 E rec (Fig. [3]). With a relative statistical un- 
certainty on the level of few percent, a Fourier analysis 
reveals both a set of frequencies of order 0(U/h) and 
0(2U/h), monotonically increasing with lattice depth. 
The measured two-particle Fock state energy E 2 is about 
10 ± 1 % lower than predicted by a simple single-orbital 
Bose-Hubbard model, which strongly indicates the rele- 
vance of multi-orbital effects in optical lattices (Fig. [3jA) . 

We compare our experimental data to an ab-initio ex- 
act diagonalization of a single-site multi-orbital system. 
For this system a contact interaction with an s-wave scat- 
tering length of 102 i2 a [35] was assumed, where ao 
denotes the Bohr radius. We find excellent agreement, 
particularly for deep lattices (Vh > 34 E rec ), when the 
diagonalization is performed on a Hilbert-space with 4 3 
single-particle orbitals, corresponding to the four lowest 
lattice bands in three dimensions (Fig. |3jV, black solid 
lines). For lower lattice depths (Vh < 34 E rec ) the re- 
sults of a calculation including 3 3 onsite orbitals appear 
even more suitable (Fig. [3j\, gray dashed lines). The 
remarkable congruence between theory and experiment 
suggests, that the influence and occupation of even higher 
lying orbitals is negligible. 

The strength of the effective multi-particle interac- 
tions has been iteratively derived from the measured 
multi-orbital energies E^f° according to U2 = E¥ , 

(7 3 = £MO _ 3f/2 and Ua = E MO _ qjj 2 _ Wz (Fig |£) 

The resulting effective three-particle interaction is at- 
tractive and the measured values agree well with the re- 
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FIG. 4: Global number statistics upon approaching the Mott insulator transition. (A) Multi-orbital quantum 
phase diffusion in a Vh = 40 E rec deep lattice was used after 3.3 ± 0.3 x 10 s 8r Rb atoms had been adiabatically loaded to 
lattices ranging from Vl — 3 to 13 E rec . The mean atom number of the individual traces differed by as little as ±1%. Although 
the coherence time at shallow lattices is significantly reduced, the visibility reliably shows dynamics down to the percent level 
(inset). (B) The corresponding Fourier power spectra reveal frequency contributions with up to six-particle Fock state energies. 
The peak positions agree with the theoretical prediction (dashed vertical lines) and are independent of Vl ■ Number-squeezing 
manifests itself in both reduced peak amplitudes and a narrowing of the spectra for increasing Vl ■ The solid lines show Gaussian 
fits to the peaks. 



suits obtained from exact diagonalization, indicating that 
direct three-particle interactions are negligible for our 
experimental parameters. The measured effective four- 
particle interaction strengths on a scale as low as 100 Hz 
are also in good agreement with theory, taking the exper- 
imental uncertainty into account. Our data shows, that 
the expansion provided by Eq. [T] quickly converges, of- 
fering the possibility to efficiently incorporate the effects 
of multi-orbital physics in refined effective single-band 
lattice Hamiltonians. 

In addition to its capability to measure coherent in- 
teratomic interaction strengths, multi-orbital quantum 
phase diffusion can be used to reveal the number statis- 
tics of many-body quantum states: The spectral weight 
of the detected frequencies carries information on the 



Fock state amplitudes c„ (Eq. pj. To demonstrate this, 
we adiabatically prepared 3D arrays of coherent (Vl — 3 
E rec ) to highly number-squeezed states (Vl — 13 E rec ) 
close to the Mott-insulator transition around 15 E rec , 
and employed quantum phase diffusion as a detection se- 
quence (see Fig. [4]A_) . From 3 to 11 E rec the time traces 
evolve from seemingly irregular oscillations into a clear 
beat signal since less frequencies contribute, as can be 
observed in the corresponding Fourier spectra (Fig.[4}3). 
This narrowing of the spectrum reflects a decreasing vari- 
ance of the atom number distribution caused both by 
number squeezing and the reduction of the average on- 
site density due to increased interactions when the lattice 
depth is raised. 

In conclusion, we have demonstrated how interaction 
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induced virtual transitions to higher lattice orbitals lead 
to effective coherent multi-particle interactions pQ. We 
have determined their strengths by monitoring long time 
traces of the collapse and revival dynamics of the matter 
wave field of a BEC. Our measurements explicitly demon- 
strate the coherence of the multi-particle interactions and 
the number state resolving capabilities of multi-orbital 
quantum phase diffusion. Multi-orbital quantum phase 
diffusion can serve as a precise experimental technique 
to determine renormalized interaction energies in lattice 
based quantum gases, which have been found to play a 
crucial role in more complex systems such as Bose-Fermi 
mixtures [26] [27j [36] . Our results also show that multi- 
orbital effects will play an important role in experiments 
working towards interaction induced dynamical creation 
of spin-squeezing or Schrodinger cat states with ultra- 
cold atoms. Furthermore, we envisage that the coherent 
multi-particle interactions demonstrated here could be 
employed to simulate effective field theories [1] and also 
help in realizing novel strongly correlated many-body 
quantum phases [22(53], e.g. with topological order [37] 
or exotic ground state properties |38j . 
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